Cosmology and two-body problem of D-branes 
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In this paper, we investigate the dynamics and the evolution of the scale factor of a probe Dp- 
brane which move in the background of source Dp-branes. Action of the probe brane is described 
by the Born-Infeld action and the interaction with the background R-R field. When the probe 
brane moves away from the source branes, it expands by power law, whose index depends on the 
dimension of the brane. If the energy density of the gauge field on the brane is subdominant, the 
expansion is decelerating irrespective of the dimension of the brane. On the other hand, when the 
probe brane is a Nambu-Goto brane, the energy density of the gauge field can be dominant, in which 
case accelerating expansion occurs for p < 4. The accelerating expansion stops when the brane has 
expanded sufficiently so that the energy density of the gauge field become subdominant. 

PACS numbers: 98.80.Cq 04.50.+h 11.25. Wx 



I. INTRODUCTION 

By the discovery of D-brane, not only string theory 
but also cosmology have been activated significantly. The 
Randall- Sundrum braneworld model 0, 0, Q is the sim- 
plest cosmological model which was induced by the idea 
of D-brane. In this model, the action of the brane is as- 
sumed to be the Nambu-Goto action. Cosmology with 
Born- Infold action has also been investigated in jj, U [|J 
and it is found that behavior of the gauge field confined to 
the brane is significantly different from that of a gauge 
field added to the Nambu-Goto brane. Interaction be- 
tween D-branes by R-R charge, which is absence in the 
Randall-Sundrum model, has been studied as a potential 
energy source that inflate the brane by many authors 
U la E EElEEiElHd- Fora review of cos- 
mology in the context of string theory, see, for example, 

m 

Since D-brane is a fundamental object in superstring 
theory, their two-body problem is also fundamental. 
Burgess et. al. studied motion of a probe brane 

in the background spacetime of source branes and found 
that there exist bound states of D6-brane and anti-D6- 
brane, which they called 'branonium'. Probe-brane dy- 
namics was also discussed in ^J|2(j. Recently cosmology 
on the probe brane was studied in the context of bounc- 
ing universe plf . 

In this paper, we investigate two-body problem and 
cosmology of D-branes. Basic approach is the same as 
[LsL l2l| but we take into account a gauge field confined 
to the probe brane, which was neglected in (l8ll2l]|. The 
motion of the brane causes the time evolution of the in- 
duced metric on it, which is seen as cosmological expan- 
sion or contraction by an observer living on the brane. 
In this sense, our pict ure is similar to that of 'mirage 
cosmology' [2I |H 1H El Thus, by following the mo- 
tion of the brane, we can also follow the evolution of the 
scale factor. We show that the gauge field on the probe 
brane, which has not been studied rigorously, can affect 
the behavior of the scale factor. 



The rest of this paper is organized as follows. In section 
ITT1 we review the p-brane solutions in supergravity as 
the background spacetime of the source D-branes. Wc 
consider the motion of a probe brane in this background 
spacetime in section IIIII and follow the evolution of the 
scale factor on the probe brane in section llVl In section 
lyland lVll we give discussion and summary, respectively. 



II. BACKGROUND SPACETIME 

We consider a system in which a probe D-brane (or 
anti- D-brane) moves within the background of N parallel 
source D-branes. In this section, we review the p-brane 
solutions in supergravity as the background spacetime of 
the source D-branes. Low-energy effective theories for 
superstring theories are given by supergravities, among 
which we consider only Type IIA and IIB here for sim- 
plicity. The effective actions include the metric, the 2- 
form potential and the scalar dilaton in the NS-NS sec- 
tor, (n — l)-form gauge potentials in the R-R sector and 
Chern-Simons terms. Here n is even for IIA and odd for 
IIB. 

To obtain a tractable system to study, we shall make a 
consistent truncation (see |'26) and references therein) of 
the action down to a simple system comprising only the 
metric Gmn, the scalar dilaton and a single (n— l)-form 
gauge potential .A[ n _i] with corresponding field strength 
-FLi . Then the background spacetime of source Dp-brane 
are determined by the following action in the Einstein 
Frame, 



S = 



D D xV^G 



R - \d M 4>0 



Al': 



2n\ [n] 



(1) 



where D = 10 and a = (5 — n)/2 is the dilaton cou- 
pling of the R-R field. Assuming asymptotic flatness and 
spherical symmetry in the transverse directions, flatness 
of the branes and an "electric" gauge field, the back- 
ground spacetime and gauge field for p < 6 are given 



2 



by. 



ds 2 = h-^-rt^jj^dx^dx" + h ip+1 ^ 8 6 mn dy m dy n 1 (2) 



,<t> = /j(3-p)/4 



Am 1 M 2 ---M p+1 = tM 1 M 2 ---M p+1 



(4) 



where = 0, 1, • • • ,p) and y m (m — 1, 2, • ■ ■ , £>— p— 1) 
are the coordinates parallel and transverse to the branes, 
respectively. We define the radial coordinate transverse 
to the brane as r 2 : 



my m y n and then, 



h(r) 



1 



p7-p" 



(5) 



Here k is an integration constant which represent the 
energy scale of the source branes: 



k = {2^f- p T{ 1 -^ l )9sll~ p N 



(0) 



where g s is the string coupling constant at infinity, l s 
is the string length scale and N is the number of the 
source branes. It should be noted that these solutions 
are reliable only for r l s . This is because supergravity 
is a good approximation of superstring theory only within 
this region, where the brane interactions are dominated 
by massless string states. 

The asymptotic behaviors of the gravitational field and 
the gauge field potential can be understood in terms of 
Gauss' law. Both behave asymptotically like ~ r~( 7 ~ p \ 
as expected from the Laplace equation, 



V 2 /M = 



d 2 8 - p d 
dr 2 r dr 



f(r) = 0. 



(7) 



Thus, the potential produced by a D6-brane is, like that 
of a point particle in ordinary 4-dimensional spacetime, 
~ r _1 . F or g lobal structures of these solutions, see, for 
example, |27|. 

On the other hand, there is no asymptotically-flat so- 
lution for p > 7. Hereafter we concentrate on p < 6 
cases. 



III. DYNAMICS OF PROBE BRANE 

In this section we consider the motion of a probe brane, 
which is assumed to be parallel to the source branes, 
in the background spacetime discussed in the previous 
section. The dynamics of the probe brane which has 
"electric" charge is determined by the Born-Infeld action 
(in the String Frame), 



Sbi = —T p 



<P+h 



det(g^ + 2l 2 F^) (8) 



and the interaction with the background gauge field 



A 



b+i] > 



Swz — ~qT p J A[ p+1 y 



(9) 



(3) Here g M „ is the induced metric on the probe brane, F t 



the U(l) gauge field strength confined to the brane and q 
is the R-R charge of the brane, which equals to ±1 for D- 
brane and anti-D-brane, respectively. Note that the field 
strength should be thermal in nature in order not to 
break the isotropy of the brane. Therefore, we interpret 
that F liv F tl " -> (F^F^) etc HI]. The induced metric 
on the brane in the String Frame is written as, 

d ~2 = e AH(D-2) ds 2 

= -h~ 1/2 (l - hv 2 )dt 2 + h- 1/2 S ij dx i dx j (10) 

where we took the static gauge, t = x° and i,j = 
1, 2, • • ■ ,p. Here we defined the velocity v of the brane 

as, 



2 dy m dy r ' 

'3 — Omn " 



dt dt 



(ii) 



Thus the motion of the brane in the dimensions trans- 
verse to the brane is described in terms of the radial 
coordinate r and the velocity v. 

Due to the spherical symmetry in the transverse direc- 
tion, the angular momenta of the brane are conserved. 
This shows that the motion is confined to the plane which 
is spanned by the initial position and momentum vectors. 
We will denote the polar coordinate in this plane by r 
and 9. Further, due to technical difficuluty, we treat the 
gauge field as a perturbation and consider the leading 
term. Then the total Lagrangian of the probe brane is, 



L 



-rah 1 



1 - h(r 2 + r 2 6 2 )(l + liF^F^) 



(12) 

where we have neglected an additive constant, m = 
T p J d p x is the "mass" of the brane, and dot denotes a 
derivative with respect to t. The independent variables 
are r, 9 and the gauge potential on the brane. The 
canonical momenta associated with these variables are, 



Pr 



dr 



Pa 



1 - h{r 2 + r 2 9 2 ) 
_ x dL_ 
89 

r 2 9 

l-h(f 2 +r 2 9 2 ) 
_ x dL 
dA~ 

4yd - h(r 2 +r 2 9 2 ) 
h 



= (l+ft,n. (13) 



-.(i + itF^n, (14) 



F 



iO 



(15) 
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where p l A is also conserved as we can see from the Euler- 
Lagrange equation. Thus the "electric field" F l ° can be 
written in terms of the other variables. On the other 
hand, the "magnetic field" F y are obtained from Bianchi 
identity, 



dnF vX + d u F Xfl + d x F^ = 0, 



as, 



C, 



const. 



(16) 



(17) 



Combining the above results, it follows that, 



F^F^ 



Ch, 



(18) 



where C is a constant which represents the energy scale 
of the gauge field. 

From Eq. I)13|l and (|14|l , we obtain the following useful 
relation: 



l 2 /r 2 



(1 + Ch) 2 + h{p 2 + l 2 /r 2 )' 



(19) 



where C = C'l^ is a dimensionless constant which repre- 
sents the energy scale of the gauge field in units of Z" 1 . 
Then the Hamiltonian can be written as, 

E = p r r + l9+p A A i -m^ 1 L 

{1 + {AD + C)h) (1 + Ch) + h{p 2 . + l 2 /r 2 ) 



(20) 



which gives the conserved energy. Here we took a gauge 
A = and D = 5ijp l A p J A /l&. Note that this agree with 
(2.22) of [3 in the limit of no gauge field, C, D -> 0. 
Hereafter, we set D = for simplicity, which means that 
there is only magnetic field. From Eq. (|20|l . we expect 
that the dynamics does not change so much even if there 
are both electric and magnetic fields. 

Following |is| , we define effective potential V e g for the 
radial motion as, 



Kff(r) 



E( Pr = 0) 
h- 1 



y/(l + Ch) 2 + hi 2 jr 2 - q . (21) 



The asymptotic behavior depends on the charge and the 
dimension of the brane. For p = 6, 



l ..-1/2 

V eS (r) -» { 1 + C - q + k{q - l)r~} 



for 



for r 





oo 
(22) 



For p = 5, 



Veff(r) 



For p < 4, 



Veff(r) 



i 

v^fe 



2kVk 



l + C-q 



^r 2 



for r 



j(T+c] + - i) 

for r 





-2 




(23) 



(24) 



2(1+C)' 



As is pointed out in |l8(, there exist stable bound orbits 
in the case of anti-6-brane. 

Behavior of the effective potential is shown in Fig. Q 
El and 01 Fig. ^ an d El shows the effective potential of p- 
brane and anti-p-brane for various p, respectively. From 
this, we can see that there can be stable bound state in 
the case of anti-6-brane, as is expected. Note that the 
position of the potential minimum, r m j n , depends on the 
angular momentum I and r m i n can be much larger than l s 
if I is sufficiently large. Fig. El shows the effective poten- 
tial of 6-brane for various C . As can be seen, qualitative 
feature does not depend on C. 

Using Eq. (T3J|, {Tl| and l|2U|). r can be expressed in 
terms of r, E,l: 



1 



r 2 (l + Ch) 2 +hl 2 



r 2 (Eh + qf 



(25) 



We can follow the motion of the brane by integrating this 
equation. Since, as can be seen from Eq. ^3 the scale 
factor on the brane is a function of r, its evolution can 
also be calculated from this equation as we discuss in the 
next subsection. Note that this equation corresponds to 
Friedmann equation and that this reduces to the Fried- 
mann equation in |2l| in the limit of C — * 0. 

The brane trajectory can be calculated as follows: de- 
fine u = 1/r and then, 



du 
d6 



i dr 

de 



Pr_ 



(26) 



Eliminating p r from Eq. I|2(J|I using this equation, we 
obtain, 



E = h- 



^(l + Ch) 2 + hl 2 {u 2 +u' 2 )-q , (27) 



from which the orbit is obtained as, 

Jl/r 

where, 



+ Bu 7 -p - u 2 ' 



A = r 2 {E 2 + 2Eq- C 2 -2(7) 
B = E 2 - C 2 . 



(28) 



(29) 
(30) 



Thus the orbit of the probe brane is equivalent to that of 
a classical nonrelativistic particle in the central potential 
proportional to r p ~ 7 , even when there exists a gauge field 
on the brane. Especially for p = 6, the bound orbit is 
closed. 
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IV. COSMOLOGY ON PROBE BRANE 

A. Evolution of Scale Factor 

From the induced metric on the brane Eq. (JTUJ, the 
scale factor a is given by, 



a = h- 1 ' i . 



(31) 



On the other hand, the cosmological time r on the brane 
is expressed as, 



= / h- l / A Jl-h{P + r 2 6 2 )dt 



h-y 



Eh + q 



h 1 ' 4 - 



1 + Ch 



=dr(32) 



y/(Eh + q) 2 - (1 + Ch) 2 - hl 2 /r 



Here we used Eq. (|19f) and (120ft in the second equation 
and Eq. H25[) in the last equation. Thus, from Eq. (|3 1|) 
and H32f> . the scale factor a can be obtained as a function 
of r. 

Here we define two characteristic radii: gravitational 
radius r g and gauge-field radius r c . The former corre- 
sponds to the Schwarzshild radius, 



It should be noted that, 

k/r 7 ~ p for r <C r g 



h(r) 



1 



for r 3> r s . 



(33) 



(34) 



The latter represents the radius, below which the approx- 
imation of the Lagrangian H12[l breaks down (l^F^F^ = 
Ch{r c ) = l): 



Ck 



V(r-p) 



i-c, 

(Cjfc)V(7-P). 



(35) 



Hereafter we consider the case r g ^> Z s because otherwise 
the scale factor does not change so much in the region 
where the background solution is reliable (r ^> l s ). 

Then let us consider a situation that the probe brane 
goes away from the neighborhood of the source branes 
(but r = ro 3> l s , r c , of cource) to infinity. When r<r g , 
the relation between the scale factor and the cosmological 
time is simple. In this case Eq. (|32[1 becomes, noting that 
h > 1 and Ch < 1, 



1/i l + Ch 



h 1 ' 4 



Eh 



dr 



E- 1 ^ 4 r^ 7 -P^ 4 (l + Ckr p - 7 )dr. (36) 



Note that this is independent on q in this limit. First we 
consider the case C — 0. Then, 



£-l fc -3/4/ r (25-3p)/4 



25 -3p 

At late time (r 3> tq), we obtain, 

„(25-3p)/4 



(25-3p)/4^ 
'0 J 



t (x r 



(37) 



(38) 



from which the evolution of the scale factor is obtained 
as. 



a(r) = h- 1 ' 4 « (kr p - 7 )-^ 4 

„ _.(7-p)/(25-3p) 



(39) 



Here (7-p)/(25-3p) = 1/7, 1/5, 3/13, 1/4, 5/19, 3/11 for 
p = 6, 5, • • • , 1. Although the expansion becomes faster 
with smaller p, acceleration phase cannot be realized. 
If C ^ 0, a correction term is added, 



o(t) cx ( t ( 7 ^)/( 25 " 3 p) - ACoi(t)), 



(40) 



where A is a constant which depends on E, k,p and a\ 
is, to leading term, 

{ r -3(7-p)/(25-3p) f or p > 4 

r- 3 / 4 logr forp = 3 

t- 1 ^ (3_p)/4 - B r -16/(3-p)(25-3p)) for p < 2j 

( 41 ) 

where is also a constant which depends on E,k,p. It 
should be noted that the effect of the gauge field decrease 
as the brane expands since its energy density decreases 
as h cx a~ 4 . 

When r becomes much larger than r g , the scale factor 
stops to evolve and becomes almost unity. The behavior 
of the scale factor in the case of no gauge field is shown 
in Fig. 0| As is expected, the scale factor evolves by 
power-law and then decelerate quickly to become unity. 
Fig. |S]shows the effect of the gauge field on the brane. As 
can be seen, the effect is very small even if C is as large 
as possible and the late-time behavior is independent on 
the existence of the gauge field. 

We can also see the evolution of the scale factor by the 
effective Friedmann equation which can be derived from 
Eq. (gSJ: 



(7-p) 2 , 2_ 2(ll-p) 

= —k 7 -pa 7 -p (1 



16 

x(l + Ca- 4 )- 2 
x [{E 2 - C 2 )a~ 4 

—Pkr^aT 7 



{\~a 4 )-^i~ 



2{Eq - C) 
-{l-a 4 )^l 



(42) 



which agrees with 22] in the limit of C — * 0. 



B. High Energy Limit 

Here we consider the probe brane to be a Nambu- 
Goto brane with gauge field and the same R-R charge 
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as a D-brane, for which the Lagrangian l|12|) is exact. In 
this case, we can take a high-energy limit (Ch 3> 1). 
Although, for a D-brane, this limit is in contradiction 
to the approximation which we used to derive (|12[1 . we 
could still obtain tendency of high-energy effect, as is of- 
ten done in higher derivative theory. In this regime, Eq. 
(|32[1 for r <C r g becomes, 



C 



For p > 4, 



Ve 

Ck 1 ' 4 
VE 2 - C 2 



4 Ck 1 ' 4 
P - 3 V^ 2 - C* 2 

r »5 oc r (p-3)/4 



(r 



-C 2 

r (p-7)/4 dr _ 
(p-3)/4 _ _>-3)/4 



(43) 



(44) 



then, 



where (7 — p)/(f> — 3) = 1/3, 1, 3 for p 
accelerating expansion is realized for p 
P = 3, 

Cfc 1 / 4 r 
T = = log — , 

VE 2 - C 2 r ' 



(45) 

6,5,4. Thus 
4. Next, for 



(46) 



then. 



a(r) 



-1/4 



r exp 



VE 2 - C 2 
Ck 1 /* ' 



(47) 



Thus the scale factor increase exponentially. Finally for 
P<2, 

T = -± SfL=(r- (3 - p)/4 - (48) 



3~PVE 2 -C 2 



then, 



a(r) = k- X ' A r - (3 - p)/4 



3 - p VE 2 - C 2 

~ Ck 1 / 4 



7-P 
" 3-p 



(49) 

It can be easily shown that the expansion is accelerating 
in this case. These analyses are confirmed in Fig. |SJ 

As stated in the previous subsection, the energy den- 
sity of the gauge field decrease as the brane expands. 
With the parametrization in Fig. H3 the gauge field is 
dominant for all over the evolution since C is sufficiently 
large so that Ch at infinity is still large (Ch(r — oo) = 
C = 1). If C is smaller than unity, late phase will behave 
like that of the case discussed in the previous subsection, 
even if accelerating expansion occurs in early phase. In 
Fig. [7| we show the cases with intermediate C. We can 
see the transition from accelerating phase to decelerat- 
ing phase. Of course, the transition occurs earlier with 
smaller C. 



C. Einstein Frame 

Finally, we give the evolution of the scale factor in the 
Einstein Frame. The procedure is almost the same as in 
the String Frame. The induced metric in the Einstein 
Frame is, 



ds 2 = -fc-(*-p)/8(i - hv 2 )dt 2 + h- ( - r - p ^ s 6 ij dx i dx j 
Then the cosmological time is, 



h-( 7 -p)/ 16 y/l- hv 2 dr. 



(50) 
(51) 



With C = 0, the scale factor evolves as, for r 3> r g , 

o(t) - h-^/ 16 cx T (T-P)*m-P)\ (52) 

where the index is (7 — p) 2 /(ll — p) 2 = 
1/25,1/9,9/49,1/4,25/81,9/25 for j> = 6, 5, ---.l. 
In high-energy limit (Ch ^> 1), for p ^ 3, 



i(r)<xt^ 2 /^ 2 , 



(53) 



where (7-p) 2 /(3-p) 2 = 1/9, 1,9, 25, 9 for p = 6,5,4,2, 1. 
For p = 3, 



/ VE 2 - C 2 \ 



(54) 



Thus, the condition that accelerating expansion occurs is 
the same as in the String Frame. It should be noted that 
the induced metric in the String Frame and the Einstein 
Frame coincide with each other for p = 3 because the 
dilaton J^J is constant in this case. 



V. DISCUSSION 

In the previous section, we dealt with a simple situa- 
tion that a probe brane goes away from the neighborhood 
of the source branes to infinity. If the probe brane ap- 
proach the source branes, the scale factor decreases as 
the inverse of that in the previous section. Then the 
other situations, for example, scattering and bound state 
of branes, are easy to imagine. In the former case, the 
brane contracts first, then bounces and finally expands. 
In the latter case, the brane continues to expand and 
contract periodically. 

In this paper, we followed the dynamics of a probe 
brane, that is, we neglected the back-reaction. This is 
justified if the probe brane is light compared to the source 
branes. This means N 3> 1, which we assumed in the 
analyses in section llVl If N ~ 1, we have to treat both 
branes equally and the self-gravity of the branes must be 
taken into account [2^, [3(j • 

Our analysis assumes a stability of the probe brane. 
There are possible instabilities due to brane bending and 
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radiation from the brane ^| . Ref. 0] has given a pre- 
liminary analysis on such instabilities. They have found, 
for p — 6, that the brane is stable classically against 
bending and that the radiation is dominated by the one 
into the bulk dilation field, which can be made sufficiently 
small by appropriate choice of the string coupling con- 
stant. One of the key assumptions they made is to treat 
the brane motion as non-relativistic one. In other words, 
their results have been obtained in the large-separation 
limit: fc/r 7 ~ p <C 1 (or r ^> r g in our notation). We 
obtain some of the expressions for the scale factor for 
r -C r g where the relativistic treatment is necessary in 
strict sence but we expect that the relativistic corrections 
do not change the stability discussed in 0|- Recently it 
was shown in |3l| that time- variations in the background 
moduli fields generally preclude the existence of stable 
elliptical orbits. 

Finally, although our study is based on the approxi- 
mated Lagrangian ll2l it would be quite interesting and 
important to study the exact Lagrangian. This will be 
our future work. 



VI. SUMMARY 

In this paper, we investigated the evolution of the scale 
factor of a probe Dp-brane which move in the background 



of source Dp-branes. When the probe brane move away 
from the source branes, it expands by power law, whose 
index depends on the dimension of the brane. If the en- 
ergy density of the gauge field on the brane is subdom- 
inant, the expansion is decelerating irrespective of the 
dimension of the brane. On the other hand, when the 
probe brane is a Nambu-Goto brane, the energy density 
of the gauge field can be dominant, in which case ac- 
celerating expansion occurs for p < 4. The accelerating 
expansion stops when the brane has expanded sufficiently 
so that the energy density of the gauge field become sub- 
dominant. Although this is not the case with a probe 
D-brane, we could still obtain tendency of high-energy 
effect of the Born-Infcld action. 

The system which is investigated in this paper is too 
simple to be our universe. However, further investigation 
will give understanding for the relation between super- 
string theory and our universe. 
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FIG. 1: Effective potential V e « for the radial motion of the probe brane, varying its spatial dimension p. Other parameters are 
set as k — I — 1 and C — 0. 




FIG. 2: Effective potential V a g for the radial motion of the probe anti-brane, varying its spatial dimension p. Other parameters 
are set as k = I — 1 and C — 0. 
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FIG. 3: Effective potential V e s for the radial motion of the probe anti-6-brane, varying the energy scale C of the gauge field 
on it. Other parameters are set as k = I = 1. 




log 10 x 

FIG. 4: Evolution of the scale factor a(r) without the gauge field on the brane for various p. Other parameters are set as, 
k = 10 s , E = 10 s , I = 10, q = -l,ro = 1. 
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FIG. 6: Evolution of the scale factor a(r) of the brane dominated by the gauge field for various p. Other parameters are set 
as, k = 10 8 , E = 10 3 , 1 = 10, q = -1, r = 1, C = 1. 
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FIG. 7: Evolution of the scale factor a(r) of the brane dominated by the gauge field for various C. Other parameters are set 
as, p = 3, k = 10 s , E = 10 3 , 1 = 10, q = -l.ro = 1. 



